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ABSTRACT 

We compute the rate of diffusive nuclear burning for hydrogen on the surface of a "magnetar" (Soft 
Gamma- Ray Repeater or Anomalous X-Ray Pulsar). We find that hydrogen at the photosphere 
will be burned on an extremely rapid timescale of hours to years, depending on composition of the 
underlying material. Improving on our previous studies, we explore the effect of a maximally thick 
"inert" helium layer, previously thought to slow down the burning rate. Since hydrogen diffuses faster 
in helium than through heavier elements, we find this helium buffer actually increases the burning 
rate for magnetars. We compute simple analytic scalings of the burning rate with temperature and 
magnetic field for a range of core temperature. We conclude that magnetar photospheres are very 
unlikely to contain hydrogen. This motivates theoretical work on heavy element atmospheres that are 
needed to measure effective temperature from the observed thermal emission and constrains models 
of AXPs that rely on magnetar cooling through thick light element envelopes. 
Subject headings: diffusion — stars: abundances, interiors, magnetic fields — stars: neutron 



1. INTRODUCTION 

The surface composition of young, isolated neutron 
stars is not well constrained by theory. Since the spec- 
trum is formed in the outermost ~ 10" 17 A/ Q , it is diffi- 
cult to predict the surface composition from supernova 
simulations, especially in light of uncertainties in the 
amount of fallback, nuclear burning in the nascent star, 
and subsequent accretion and spallation reactions. On 
the observational side, the absence of well-understood 
atomic spectral lines, and the extreme physical condi- 
tions at the photosphere impede progress using spec- 
tral modeling, though some inferences have been made ( 
Pavlov, Sanwal k, Zavlin 2002). 

Diffusive nuclear burning (DNB) is a mechanism by 
which the composition of young neutron star (NS) sur- 
face layers can evolve with time (Chang & Bildsten 2003; 
2004, and Chang, Arras & Bildsten 2004, hereafter pa- 
per I, II and III respectively). The simple picture of 
DNB (Chiu & Salpeter 1964, Rosen 1968) can be under- 
stood as follows (see Fig. 1 of paper I). Rapid sedimenta- 
tion (~ 1 sec) allows H to settle above heavier elements. 
While nuclear burning timescales for H at the surface 
are far too long to be of interest, a diffusive tail readily 
extends to depths at which proton captures occur, in- 
creasing the burning rate by many orders of magnitude. 

Anomalous X-Ray Pulsars (AXP's) and Soft Gamma- 
Ray Repeaters (see e.g. Woods & Thompson 2004 for 
a review), are a class of strongly magnetized (B ~ 
10 14 " 15 G), hot (T bb - (5 - 7) x 10 6 K) 3 , young 
(P/2P ~ 10 3 ~ 4 yr), and slowly rotating (P ~ 5 - 10 sec) 
neutron stars. We will refer to them as "magnetars" , 
following Duncan and Thompson (1992). In this Let- 
ter, we show that the high temperatures and substan- 
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tial alteration of the equation of state (EOS) by super- 
strong magnetic fields combine to increase the rate of 
DNB by many orders of magnitude as compared to other 
observed neutron stars. This conclusion is robust for a 
variety of heavy proton capturing elements, even when a 
He layer exists between H and the proton-capturing el- 
ement. Consequently, magnetar photospheres are most 
likely composed of heavy elements, possibly helium. 

2. REVIEW OF DNB 

The rate at which a column yu — J pndz of H is 
burned by the proton-capturing substrate is (Papers I 
and II) 

Vh = = dz — — , (1) 

Tcoi J TH{nH,n pc , 1 ) 

where z is the depth, uh and n pc are the number den- 
sity of H and proton-capturing substrate, respectively, 
T is the temperature, and th — {(&v)n pc )~ 1 is the lo- 
cal lifetime of a proton. The two key ingredients in this 
formula are the temperature profile and concentration, 
fn = nn/n-pc, of H in the diffusive tail. The temperature 
profile for strongly magnetized NS envelopes is reviewed 
in Ventura & Potekhin (2002). 

The concentration of a trace particle with charge Z t 
and mass number A t in a background Zb and Aj, de- 
pends on the thermal (kgT), Fermi (Ep), and Coulomb 
interaction (£c) energies. In the non-degenerate limit, 
k B T ^> Ep, the energy gained by separating the two 
species across a pressure scale height is ~ ksT. The 
resulting concentration profile / = nt/nb is a power- 
law with column (paper I) / cx y s with exponent 8 = 
At(Zb + 1)/Ab — Zt — 1. In the degenerate limit, if 
Z t /A t ^ Zb/Ab, the separation energy is ~ Ep over a 
scale height, leading to an exponential concentration pro- 
file (paper I) / oc exp \-{Z t - Z b A t /Ab)E F /k B T}. How- 
ever, if ksT -C Ep and Z t /A t — Zb/Ab, the dominant 
separation energy is provided by Coulomb interactions 
Eq ~ Z 5 / 3 e 2 /et e where the mean electron spacing is 
a e = (47rn e /3) _1 / 3 . The resulting concentration profile 
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is (paper III) / oc exp -0.9Z t (^ /3 - Z t 2/3 )e 2 /a e k B T 

The composition profile drops off much more rapidly in 
the degenerate limit (exponential) as compared to the 
nondegenerate limit (power-law). 

For hot NS's (T e > 10 62 K) such as magnetars, DNB 
is "diffusion limited" (see paper II), meaning the burn- 
ing occurs in a thin layer at which the nuclear burning 
time becomes comparable to the downward diffusion time 
th = Tdifi breaking the assumption of diffusive equilib- 
rium (paper I, II). The hydrogen concentration profile is 
rapidly cut off in this region. Cool NS (T e < 10 62 K) 
have "nuclear limited" DNB with th < Tdiff, even in the 
burning layer. In this case, burning occurs in a degener- 
ate layer where the temperature profile becomes isother- 
mal and hydrogen fraction is exponentially dropping (pa- 
per I). 

Paper II studied the effect of a "standard" pulsar 
strength (B — 10 12 G) magnetic field on the composi- 
tional structure. The effect of electron quantization into 
Landau levels was taken into account in the opacity and 
EOS. For B = 10 12 G, it was found that the burning rate 
decreased by more than an order of magnitude. Most of 
this decrease is due to the smaller temperature at fixed 
column for a given core temperature T c . In addition, fn 
in the burning layer decreased as compared to the B = 
case. The reason, clearly seen in Fig. 4 of Paper II, is 
that the electric field above the burning layer is on av- 
erage larger for B = 10 12 G than for B = 0, leading to 
a smaller diffusive tail of H. However, this is due to an 
unlucky coincidence of the lowest Landau level nearly co- 
inciding with the burning layer, so that the electric field 
shows oscillations due to the De Haas- van Alphen effect. 
Larger magnetic fields (see the B = 10 13 G line) cause 
the electric field above the burning layer to be weaker 
than the B = case, increasing yn in the burning layer. 
As a result, DNB can be much faster in magnetars than 
for standard radio pulsar fields. 

3. MAGNETIC FIELD DEPENDENCE OF DNB 

Magnetic fields significantly modify the opacity and 
EOS when electrons occupy only the ground Landau 
level, i.e. fuo c ^ e > max(fcsT, Ep) where huj c ^ e = 
ll.57.Bi2 keV is the cyclotron energy. There are two main 
effects. First, electron scattering and free-free Rosse- 
land opacities are reduced by a factor ~ (ksT [ /ftw Cye ) 2 . 
Hence, for fixed T c , larger B implies lower temperature 
at a fixed column (paper II, Potekhin et al. 2003), 
and a smaller burning rate. The photosphere also 
moves to higher density, and the T e (T c ) relation be- 
comes highly dependent on B and 9, the angle between 
the field and normal. Second, the gas becomes more 
nondegenerate. The ratio of the B = Fermi energy 
Ef,3D to the ground Landau level Fermi energy Ep,iD 
is Efsd/Ef^d ~ (E F , 1D /tiuJc,e) 2/3 < 1, pushing the 
point ksT = Ep deeper into the envelope. This dramat- 
ically alters the composition profile, allowing the slower 
power-law decrease all the way to the burning layer, in- 
creasing the burning rate. In effect, the parameter space 
for "diffusion-limited" DNB is made larger. 

The opacity change decreases burning while the degen- 
eracy change increases burning. Which one wins? The 
answer is that in the high temperature limit, the burning 
rate decreases with B, and in the low temperature limit 
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Fig. 1. — Hydrogen column lifetimes for p — 12 C as a function of 
T c for envelopes with vertical magnetic field of different strength 
and yn = 10 4 gcm -2 . We show the low field case (B = 10 9 G, 
thick-solid-line), pulsar field strength (B = 10 12 G, dotted-line, 
and B = 10 13 G, short-dashed-line), and magnetar field strength 
represented by the long-dashed line (B = 10 14 G). The bottom 
panel shows T e as a function of T c for the different field strengths. 



it increases with B (in the large B limit), as we show in 
this section. 

We extend the methodology and results of paper II 
to higher magnetic fields and core temperatures. We 
calculate DNB timescales assuming a vertical magnetic 
field, 12 C as the proton-capturing element, a hydrogen 
column of yn = 10 4 g cm~ 2 , and a fixed surface gravity 
g = 2.43 x 10 14 cm s~ 2 . Varying the magnetic field di- 
rection was discussed in Paper II. The ions are assumed 
to be an unmagnetized ideal gas. We use the results 
of Chabrier & Potekhin (1998, 2000) for the electron 
equation of state and Potekhin et al. (1999) for the 
electron conductivity. We ignore Coulomb interactions 
in this section for simplicity, but include it in the next 
(see Paper III for details). Radiative opacities are from 
Potekhin & Yakovlev (2001). The composition profile is 
computed including diffusion and nuclear burning (see 
paper II). The ion-ion diffusion coefficient in the nonde- 
generate regime is from Alcock & Illarionov (1980) and in 
the liquid regime from Brown, Bildsten & Chang (2002). 
We refer the interested reader to paper I and II for ad- 
ditional details. 

Fig. ^ shows the DNB timescale, r co i = yn/yn as a 
function of T c for a range of B at fixed yn ■ The low field 
(B = 10 9 G) case (solid line) is included as a baseline. 
In the low T c limit, t co i initially increases with B until 
B ~ 10 12 ~ 13 G, after which it decreases with B. In the 
high T c limit, t co \ increases with B. The dependence on 
core temperature is power law for high T c and exponen- 
tial for low T c . The rate at which hydrogen is depleted 
from the photosphere, y p h ~ 500 g cm _2 r e 5 g 4 Bi 5 , can 

5/12 

be found using the scaling relation t co \ oc y H . The 
column lifetime at the photoshere is far shorter than the 
magnetar age (~ 10 3 ~ 4 yr). 
We estimate the burning rate in the high temperature 
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limit following paper II. Using free-free opacity with elec- 
trons in the ground Landau level and ideal gas pressure 
from electrons and ions gives the temperature profile 
(Ventura & Potekhin 2001) T 7 = l.l(T e 2 6 y 6 /Bi 5 ) 4 / 13 . 
The diffusion coefficient for nondegenerate ions gives a 
diffusion time T diff ~ 2.5 x I0 4 sec( J Bi 5 /T e 2 6 )T 7 9/4 . The 
lifetime of a proton to nonresonant capture on C is 
t h ~ 1.4 x 10 15 sec(T e 2 6 /Bi 5 )T 7 - 23 , expanded around 
T-j = 1. Setting r^g = th gives the temperature, concen- 
tration, column, and capture rate in the burning layer. 
Placing these results into eq. Q gives the rough analytic 
formula 

5/12 / „ \ 1.2 

t co i ~ O.lyr 



10 4 g cm 2 



Vh 



D 



10 15 G 



10 6 K 



2.4 

(2) 



in the high T c limit. 

In the low temperature nuclear-limited regime, we es- 
timate the burning rate following the analytic calcula- 
tion in Paper I. In this case, the burning occurs in a de- 
generate and nearly isothermal region, so the hydrogen 
fraction is f H ~ (y H /yd cg ) 17/12 exp(l - E F /2k B T c ) = 

(z/.ff/2/deg) 17/12 exp(l - (y/?/dog) 2/3 ), where the column at 
which electrons in the ground Landau become degener- 
ate is y dcg ~ 3.9 x 10 8 g c m - 2 Bi 5 r c 3/ 8 2 . The first fac- 
tor in fn accounts for the decrease in the outer non- 
degenerate envelope, while the second factor gives the 
exponential decrease in the degenerate regime. The 
temperature is nearly T c , with a small correction scal- 
ing as y _1/3 . Inserting this into the p- 12 C rate gives 
exp(-137/T 6 1/3 ) ~ exp(-137/T t 1 / 6 3 )exp[-( 2 / nuc / 2 /) 1 /3], 
where the burning rate changes over a column y nuc ~ 
7.2 x 10 7 g cm~ 2 B 15 T^ 2 6 T- 7 . Using the method of steep- 
est descents (Paper I) for the integral in eq^ the domi- 
nant scalings for the column lifetime in the low T c limit 
are 



r col oc B 15 7/12 y ff 5 / 12 exp(29.5/T f 



l/3x 
c,8 I 



(3) 

The scaling of r co i with B can be explained as two factors 
of B from the width of the burning layer and the density 
in the burning layer, and the decrease of composition in 

the nondegenerate envelope 2/ de g 7/12 oc _B 15 17/12 . 

4. HEAVIER P-CAPTURING ELEMENTS AND A HELIUM 
BUFFER 

In section|21we found that the (diffusion-limited) burn- 
ing rate for p — 12 C is fast enough to deplete H from 
the photosphere of a magnetars in a few hours. We 
now repeat the calculation for heavier elements Mg, 
Si and Fe. The Mg and Si rates are f rom NACRE 
(http://pntpm.ulb.ac.be/Nacre/nacre.htm) and Fe from 
Schatz (2004, private communication). In spite of the 
large Coulomb barrier for these elements, we find that 
H is easily consumed. For Mg, Si and Fe, the column 
lifetimes for yn = 10 4 gcm -2 at B = 10 14 G are, respec- 
tively, r col = 0.005, 0.05 and 2 years for T e = 7 x 10 6 K, 
0.01, 0.1 and 5 years for 5 x 10 6 K, and 0.03, 0.7, and 
100 years for 3 x 10 6 K. These timescales are all much 
shorter than magnetar spindown ages ~ 10 4 yrs. 

For sufficiently large Z, the reaction rate is so small 
that DNB is nuclear-limited instead of diffusion-limited. 
Column lifetimes then increase exponentially with tem- 
perature, as opposed to the slower power-law increase 
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Fig. 2. — Total H DNB rate including a He buffer between 
the H and C layers. In the lower panel, the solid lines denote 
number fraction of the various species and the dashed line is the 
temperature profile. For T e = 4 X 10 6 K and yn = 2 X 10 4 gem -2 , 
the column lifetime is r co i = 5 hours. The total He column is 
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10 8 gcm~ 2 and the B field is B = 10 14 G. 



in the diffusion limited regime. The transition between 
nuclear and diffusion limited burning is given by the 
condition r^ifF = th evaluated at the burning layer. 
Since the nuclear reaction rate is oc exp(i? /fcsr) where 
E oc Z 2 / 3 T 2 / 3 is the Gamow energy (Clayton 1983), 
and T e oc T c 1/2 (Ventura & Potekhin 2002), the transition 

temperature T e>tt oc T c 1/2 oc Z. For T e > 5 x 10 6 K, DNB 
occurs only in the diffusion limited, and hence rapid, 
regime for elements up to and including Fe. 

We now consider the effect of a He buffer between H 
and the underlying p-capturing elements. The H concen- 
tration decreases into the He layer as /# ~ (y/yH)~ 17 ^ 12 , 
and the C concentration as fc — (y/j/He) 2 in a nondegen- 
erate layer and exponentially in a degenerate layer, due 
to Coulomb interactions. Hence an arbitrarily thick He 
buffer would strongly suppress DNB. However, thermal 
stability for the 3a reaction sets a maximum He column 
2/Hc < 10 8 gcm" 2 for T > 2 x 10 s K. In Fig. |2 we plot 
the H/He/C profile for a magnetar (B = 10 14 G) with 
an effective temperature of T e — 4 x 10 6 K. Because of 
the high temperatures, even a maximally thick He buffer 
makes little difference in the column lifetime, which is 
set by the diffusion time down to the burning layer. A 
pure H/C envelope at this T e has r co i ~ 30 hours. In- 
cluding the He buffer actually reduces the lifetime to 
t co i ~ 5 hours! The higher burning rate is the result of 
faster diffusion of H through He rather than C in the dif- 
fusion limited regime. To get the rough factor, consider 
the diffusion time r^iff oc Z^A b 1/2 , where the scaling is 
found from the diffusion coefficient in the nondegenerate 
limit. We find Tdiff is shorter by a factor of 5 for He 
compared to C. Hence for the high T c in magnetars, our 
estimates for the DNB burning rate are not slowed by 
the presence of a He buffer and in the case of diffusion- 
limited DNB may be enhanced. 
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5. SUMMARY AND CONCLUSIONS 

We have shown that diffusive nuclear burning re- 
moves hydrogen from the photospheres of magnetars on 
a timescale much shorter than the age. This conclusion 
holds for a range of proton capturing elements up to Fe. 
Our new understanding of stratification in Coulomb liq- 
uids allows us to include a maximally thick He buffer, 
which increases the burning rate in the diffusion limited 
regime. 

We conclude that magnetar photospheres will not con- 
tain hydrogen unless it can be supplied at a rate 

The burning time is so short that the supply of hydro- 
gen must be continuous over the timescale of the ob- 
servations. The relativistic outflow expected from mag- 
netars likely precludes accretion of hydrogen. However, 
spallation reactions due to magnetospheric currents may 
generate fresh hydrogen. 

The absence of hydrogen will modify our understand- 
ing of both the thermal evolution and the thermal spec- 
trum of magnetars. Atomic physics (e.g. Lai 2001) 
and spectral modelling of heavy element atmospheres in 
B ~ 10 14 G fields is relatively unexplored. The light 
atmosphere models of Heyl & Hcrnquist (1997), Ozcl 
(2001) and Ho & Lai (2001) assume fully ionized H and 
may be significantly different if the photospheric material 
is He or heavy elements. Even after T e has been deter- 
mined, the question remains as to whether the data re- 
quire an additional heat source (e.g. field decay; Heyl & 
Kulkarni 1998) or whether the apparent "hotness" is due 



to a thick layer of low opacity material (i.e. helium; Heyl 
k Hernquist 1997). Models of AXP's that rely on mag- 
netar cooling with a H envelope thicker than even a pho- 
tosphere are difficult to realize given the rapid timescales 
we have found for H consumption. Thermal stability re- 
stricts hot (T > 2 x 10 8 K) He envelopes to a maximum 
column of y Ho < 10 8 gcm~ 2 . Since T c » 2 - 3 x 10 8 K 
(Yakovlev, private communication) for a non-supcrfluid 
NS between the ages of 10 3 — 10 4 years, the effective tem- 
perature of such a maximal envelope (j/hc = 10 8 g cm" 2 
with C below) is T e w 3 x 10 6 K at B = 10 15 G. In com- 
parison a thermally unstable pure H or pure He envelope 
would have T e « 5 x 10 6 K and T e « 4 x 10 6 K, respec- 
tively. Therefore, passively cooling neutron stars without 
internal heat sources (such as magnetic field decay; Heyl 
& Kulkarni 1998 and references therein; also see Arras, 
Cumming & Thompson 2004) fall short of the observed 
luminosity of AXPs by a factor of 3 to 30 (Perna et al. 
2001). 
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